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Abstract 

The pseudoperturbative shifted - I expansion technique PSLET 
[12,16] is generalized for states with arbitrary number of nodal ze- 
ros. Interdimensional degeneracies, emerging from the isomorphism 
between angular momentum and dimensionality of the central force 
Schrodinger equation, are used to construct part of the D - dimen- 
sional anharmonic oscillator bound - state spectra. PSLET results 
are found to compare excellently with those from a series [7], exact 
and an open perturbation [9] solutions. 
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1 Introduction 



Besides its importance in quantum field theory and chemical physics, quantal 
isotropic anharmonic oscillator 

V(r) = a q 2 + aq 4 , (1) 

is also a popular theoretical laboratory for examining the validity of approx- 
imation techniques [1-12]. Therefore, it has been a subject of an enormous 
number of papers ( long list of these could be found in [3]). However, most of 
these studies are devoted to one spatial dimension (ID). Nevertheless, there 
has been a considerable interest, since the end of the 1960s as a result of 
studies in the far infrared and microwave regions, in the analysis of two - 
and three - dimensional ( 2D and 3D, respectively) anharmonic oscillators 
[6, and references therein]. Among the most recent ones, to the best of our 
knowledge, exist the work of Lakshmanan et al. [6] and Taseli [7] . They have 
treated the 2D - and 3D - anharmonic oscillators via a phase - integral [6] 
(PIM) and a series solution, weighted by an appropriate function, [7] (SSM) 
methods. Hence, further studies of such oscillators should be of great interest. 

In many problems, on the other hand, the Hamiltonian does not con- 
tain any physical parameter suitable for a perturbation expansion treatment. 
More often, the Hamiiltonian contains physical parameters, but, typically, ze- 
roth - order solutions for especial values of these are not tractable or good 
starting approximations. Moreover, it is well known that the implementation 
of Rayleigh - Schrodinger perturbation theory, or even naive perturbation se- 
ries, expresses the eigenvalue of the AHO as a formal power series in a which 
is quite often divergent, or at best asymptotic, for every a / 0. One has 
therefore to resort to summation tools to sum up such series [12-15]. Hence, 
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apparently artificial perturbation theories have been devised and proven to 
be ways to make progress [1,9,10,12,16-21, etc.]. Yet, in the simplest case, 
analytical calculations can aid numerical studies in areas where numerical 
techniques might not be controlled. For example, when bound - state wave 
functions with arbitrary nodal zeros are in demand for certain singular po- 
tentials ( a next level of complexity), analytical solutions can supply a basis 
for numerical calculations. 

In a preceding paper [12], we have treated the nodeless states of the AHO 
(1) via a pseudoperturbation shifted - / expansion technique (PSLET), where 
/ denotes the angular momentum quantum number. Successfully, the same 
recipe has been applied to quasi - relativistic harmonic oscillator [16], and 
spiked harmonic oscillator, etc [17]. 

Encouraged by its satisfactory performance in handling nodeless states, 
we feel tempted to generalize PSLET recipe ( in section 2) for states with 
arbitrary number of nodal zeros, k > 0. Moreover, in the underlying time - 
independent radial Schrodinger equation, in h — m — 1 units, 



the isomorphism between orbital angular momentum I and dimensionality 
D invites interdimensional degeneracies [7,22-24]. Which, in effect, allows us 
to generate the ladder of excited states for any given k and nonzero / from 



2~cUf 



+ 



1(1 + 1) 
2q 2 



+ V(q) V k ,i(q) = E k ^ k!l (q), 



(2) 



the /=0 result, with that k, by the transcription D 
E k j(D) is the eigenvalue in D - dimensions, then 



D + 21. That is, if 



E k>l (2) = £ fc) ,_i(4) 



E k>1 (2l) = E k>0 {2l + 2) 



(3) 



for even D, and 
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(4) 



for odd D. For more details the reader may refer to ref.s [7,22,25]. We 
therefore calculate, in section 3, the energies for 2D - and 3D - anharmonic 
oscillators, for a given number of nodes k and different values of I, and 
construct part of its D - dimensional bound - state spectra. We compare our 
results with those from PIM [6], SSM [7], exact and an open perturbation 
solutions [9]. Section 4 is devoted for concluding remarks. 

2 The generalized PSLET 

Although some of the following expressions have appeared in previous articles 
[12,16,17], we would like to repeat them to make this article self contained. 

We simply start with shifting the angular momentum quantum number 
/ in (2) through I — I — (5 and use 1/7 as a pseudoperturbation expansion 
parameter. Where (3 is a suitable shift introduced to remove the poles that 
would emerge, at lowest orbital states with I — 0, in our expansions below. 
Hence, equation (2) reads 

where Q is a constant that scales the potential V(q) at large - Id limit ( the 
pseudoclassical limit [22]) and is set, for any specific choice of l D and k, equal 
to l 2 at the end of the calculations. Here l D — I + (D — 3)/2, to incorporate 
the interdimensional degeneracies associated with the isomorphism between 
angular momentum and dimensionality D. Hence, / — >■ I = Id — (3 through 
out this paper. Next, we shift the origin of the coordinate system through 
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x = l l t 2 (q — q )/Qo, where q Q is currently an arbitrary point to be determined 
below. Expansions about this point, x = (i.e. q = q ), localize the problem 
at q Q and the derivatives, in effect, contain information not only at q Q but 
also at any point on g-axis, in accordance with Taylor's theorem. Equation 
(5) thus becomes 



Id 2 a 2 ~ 



Ho 



tt fc ,,(a;) = jE k ^ Kl {x), 



(6) 



with 



i 



V(x(q)) = q 2 J 



1 V(qo) 
2q 2 Q 



+ l 1 l 2 B l x + Y j v {n \x)l 



n/2 



n=0 



(7) 



where 



v^\x) = B 2 x 2 + 2 -^.. 



(8) 



v {1 \x) = -{2/3 + l)x + B 3 x 3 , 



(9) 



v {n \x) = B n+2 x n+2 + (-1)" (2(3 + 1) + x n 

+ (_i)«^±I) ( n _i) x («-2) ; n >2, (10) 



(ii) 
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It is then convenient to expand E k ,i as 



E k , = E 4? r n 



(12) 



n=-2 



Equation (6), along with (7)-(ll), is evidently the one - dimensional Schrodinger 
equation for a perturbed harmonic oscillator 



X k (x) = X k X k (x 



(13) 



where w 2 = 2B 2 , 



So = I 



1 Q 2 V(Qq) 

2 Q 



2/3+1 (3((3 + l) 
2 + 2Z ' 



(14) 



and P(x) represents the remaining terms in Eq.(6) as infinite power series 
perturbations to the harmonic oscillator. One would then imply that 



1 

+ I 



1 « 

2 Q 

W+l) + A (o) 



2(3 + 1 n 1. 

oo 
n=2 



(15) 



and 



a, = al E 4?r (n+1) - 



n=-2 



(16) 



Hence, equations (15) and (16) yield 
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F (-2)_ 1 V(q ) 



(17) 



(-D = £ 



2/5 + 1 1. 



(18) 



W+l) + A (o) 



(19) 



4? = A^o 2 ; n > 1. (20) 



f— 2) 

Where g G is chosen to minimize E^. t , i. e. 

, '- = and — > 0. 21 

dq Q dq 2 Q 

Hereby, V(q) is assumed to be well behaved so that E\r^ has a minimum q Q 
and there are well - defined bound - states. Equation (21) in turn gives, with 

i = VQ, 



lD-P = y/(gV'(qo). (22) 

Consequently, the second term in Eq.(7) vanishes and the first term adds a 
constant to the energy eigenvalues. It should be noted that the energy term 
PEf~^ corresponds roughly to the energy of a classical particle with angular 
momentum L z —l executing circular motion of radius q Q in the potential V(q a ). 
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It thus identifies the zeroth - order approximation, to all eigenvalues, as a 
classical approximation and the higher - order corrections as quantum fluc- 
tuations around the minimum q , organized in inverse powers of I. The next 
correction to the energy series, lE^] 1 ^, consists of a constant term and the 
exact eigenvalues of the harmonic oscillator w 2 x 2 /2.The shifting parameter 
(3 is determined by choosing lE^^—O. This choice is physically motivated. 
In addition to its vital role in removing the singularity at I — 0, it also re- 
quires the agreements between PSLET eigenvalues and eigenf unctions with 
the exact well known ones for the harmonic oscillator and Coulomb poten- 
tials. Hence 



P = - 



2 + {k + jw 



w 



1 



qoV"(qo) 

V'(q ) 



(23) 



where primes of V(q ) denote derivatives with respect to q Q . Then equation 
(6) reduces to 



1 rl 2 00 



E 

Ln=l 



2 j?{n-i)j- 

%E kl I 



(24) 



Setting the wave functions with any number of nodes as 



Vk,i(x(q)) = Fk,i( x ) exp(U k j(x)), 



(25) 



equation (24) readily transforms into the following Riccati equation: 



F k ,i(x) 



1 oo 

— 9 (U'^x) + U' Kl {x)U' Kl {x)) + £ v^(x)l-^ 



n=0 



E^r^ n 



F' k ^x)U' Kl {x) - ^» = 0, 



(26) 
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where the primes denote derivatives with respect to x. It is evident that this 
equation admits solution of the form 



7t(")/ 



n=0 



n=0 



(27) 



oo k— 1 



n=0p=0 



(28) 



where 



rt+l 

E 

m=0 



/-) /rra-l . n — n 

J -'m,n,k i J - / 0,n,fc u ; 



(29) 



n+l 

4 n) (*) = E 

m=0 



2m 



(30) 



Substituting equations (27) - (30) into equation (26) implies 



F k ,i(x) 



- E ( u l n) J ~ n/2 + G< k n) r (n+1)/2 ) 



n=0 



-i oo n 

± \^ (jj( m )jj( n - m )j-n/2 , Q(m)Q(n-m)j-{n+2)/2 

c\ / j / j \ k k k k 



n=0 m=0 



oo oo 

+ 2^ m) Gi n - m) z"-( n+1 )/ 2 ) + jr «( n »r/ 2 - jr i 2 o4T 1)J " n 



n=0 



n=l 



g^(n) / --„/ 2 + G (n) [ -(„ + l)/2^ 



.71=0 



" ^» = (31) 



The above procedure obviously reduces to the one described by Mustafa and 
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Odeh [12,16,17], when k = 0. Moreover, the solution of equation (31) follows 
from the uniqueness of power series representation. Therefore, for a given k 
we equate the coefficients of the same powers of I and x, respectively. For 
example, when k — 1 one obtains 

D lfill = -w, Ui°\x) = -wx, (32) 

^1,0,1 — , o 01 — , [66 

w ' w 

Co.o.1 = ^ (2C 1)0l i + 2/3 + 1), (34) 

^ = ^(%-^). <»> 
A A i = ^ (J + C ,o,i C 1|0 ,i - \{2(3 + 1)) , (36) 

^ " g ^ 2 + 0,1 2 2~J ' (37) 

etc. Here, we reported the nonzero coefficients only. One can then calcu- 
late the energy eigenvalues and eigenf unctions from the knowledge of C mjTlt k, 
D m ,n,k, and affi. in a hierarchical manner. Nevertheless, the procedure just 
described is suitable for a software package such as MAPLE to determine 
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the energy eigenvalue and eigenfunction corrections up to any order of the 
pseudoperturbation series (12). 

Although the energy series, equation (12), could appear divergent, or, at 
best, asymptotic for small I, one can still calculate the eigenenergies to a very 
good accuracy by forming the sophisticated [N,M] Pade approximation [1] 

P#(l/Z) = (P + PJI + ■■■ + Pm/T M )/{1 + qi/l + ■■■ + q N /l N ) 

to the energy series (12). The energy series (12) is calculated up to Ef\ /I s 
by 

Ek,i = + Ef} + ■■■ + E^/P + 0(1/P), (38) 



and with the Pf{l/l) Pade approximant it becomes 

E k MA]=PEif + P 4 \l/l). (39) 



Our recipe is therefore well prescribed. 



3 D - anharmonic oscillator spectra 

In this section we consider the phenomenologically useful and methodically 
challenging AHO interactions (1). and illustrate the above mentioned pro- 
cedure. 

The substitution of (1) in (23), for k > 0, implies 



w 



8a q + 2Aaql 
\ 2a q + Aaql ' 



(40) 
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and Eq.(22) yields 



Id + \ I l + (2fc+l)^ 



8a q + 2Aaq^ \ 2 r — — - 

2 aoqo + A aq l ) =q °^ a ° + i °">°- (41) 



Once q Q is determined (often numerically) the coefficients C mtn ^, E mn ^ and 
a^l are obtained in a sequential manner. Then, the eigenvalues, equation 
(38), and eigenfunctions (25), along with (27)-(30), are calculated in one 
batch for each value of k, D, I, a a , and a. 

Tables 1 and 2 show the 2D - and 3D - AHO energies, respectively, for 
different values of k, I, and g = a/2. PSLET results, Ep, compare excel- 
lently with those reported by Taseli [7] via a series solution weighted by an 
appropriate function, Ess- The 3D - AHO energies for k — 1 and / = 0, or 
equivalently ID - AHO third excited state, are displayed along with those 
of Bessis and Bessis [9], E BB , via an open perturbation technique and the 
exact ones, E ex , using Bargman representation [26] ( direct numerical integra- 
tions) for different anharmonicities. Obviously, our results are in quantitative 
and/or qualitative agreements with the other ones. They are also in good 
agreement with those reported by Kleinert [3] via a variational perturbation 
approach. 

Clearly, the accuracy of PSLET increases for larger k and/or I. The 
P£(l/l) Pade approximant enhances the accuracy, although its effect is not 
dramatic for weak anharmonicities. Hence, one proceeds with confidence 
and obtain, via (39), the 3D - energies ( table 4) for states with k — 1 and 
I — 1, 5, 10. In doing so, one should keep the stability of the Pade squence in 
point ( for more details on this issue the reader may refer to [16]). The same 
recipe is used to calculate the 2D - AHO energies ( table 5) for states with 
k — 0, 1 and I = 0, 1, 5, 10. Then the D - dimensional spectra follow from the 
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implicated wisdom of (3) and (4) (i.e.: the 2D - and 3D - AHO energies are 
the basic ingredients for the construction of the D - dimensional ones [22]). 

4 Concluding remarks 

In this work we have generalized PSLET [12,16,17] for states with arbitrary 
number of nodal zeros, k > 0. Starting with the "radical" central force 
problem, represented by the radial Schrodinger equation, and generalizing 
the angular momentum to the D - Dimensional one ( i.e.: I — > Id = l + (D — 
3)/2), we have treated the AHO in D - dimensions. The comparison between 
PSLET results with the other ones, including direct numerical integration, 
is readily satisfactory. 

Although we have used Pade approximants to improve the numerical per- 
formance of PSLET, it is by no means clear whether the Pade approximants 
are necessarily the most effective ones. It has been suggested by Weniger [14] 
that much better results could be obtained via a class of sequence transfor- 
mations [27] ( The details of which could be found in [27-30, and references 
therein]). However, such interesting investigations already lie beyond the 
scope of the attendant proposal. 

Finally, the applicability of PSLET extends far beyond the present D 
- dimensional anharmonic oscillator model. It could be applied to angular 
momentum states of multi - electron atoms [23,24], quark - antiquark models 
[31], 2D - Hydrogenic donor states in an arbitrary magnetic fields [32], 2D - 
magnetoexcitons [33], two - electron quantum dots [34], etc. 
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Table 1: 2D anharmonic oscillator (1) energies, in h — m — 1 units. Where 
a = 1/2 and g = a/2 to recover Taseli's results, E S sm, via series solution 
[7]. E P represents PSLET results, Eq.(38), and £[4,4] shows the P 4 4 (l/Z) 
Pade approximant,Eq.(39). 



9 


k 


I 


2£ P 


2£[4,4] 


Essm 


1(T 4 








2.000199955022 


2.000199955022 


2.000199955022 






5 


12.00419695953 


12.00419695953 


12.00419695953 




1 


1 


8.002398591662 


8.002398591662 


8.002398591662 




3 




16.00958904459 


16.00958904459 


16.00958904459 


1 








2.947835 


2.952052 


2.952050 






2 


10.390626203 


10.390627276 


10.390627295 






4 


19.217523488 


19.2175234955 


19.2175234959 




1 


1 


15.48277174 


15.48277148 


15.48277158 


10 4 








50.75164 


50.54788 


50.54804 






4 


368.030083 


368.030082436 


368.030082448 




1 





205.3783 


205.3774 


205.3777 






2 


394.577414 


394.577403 


394.577407 
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Table 2: Same as Table 1 for 3D anharmonic oscillator (1) energies. 



9 


k 


I 


2E P 


2£[4,4] 


Essm 


1(T 4 








3.0003748969 


3.0003748969 


3.0003748969 






10 


23.014356719 


23.014356719 


23.014356719 


1 







4.648511 


4.648815 


4.648813 






1 


8.380337 


8.38034245 


8.38034253 






5 


26.528917558 


26.528917558 


26.528917558 




1 


3 


27.898417763 


27.898417756 


27.898417760 


20 





1 


19.783266 


19.7832518 


19.7832519 






2 


30.057200 


30.057199029 


30.057199045 






5 


65.961500037 


65.96150003049 


65.96150003068 




1 


1 


44.209282 


44.209279007 


44.209279973 


10 3 








38.092 


38.086822 


38.086833 







3 


149.439046 


149.439045568 


149.439045581 
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Table 3: Three-dimensional energies for the k — 1 and I = state or 
equivalently one- dimensional third excited state energies for V(q) = y + ag 4 . 
Ebb denotes Bessis and Bessis results [9] and the exact ones E exact , reported 
therein, from direct numerical integrations for different anharmonicities. 



a 


Ep 


£[4,4] 


E 


Ebb 


0.002 


3.53674413 


3.536744133 


3.53674413 


3.53674 


0.01 


3.67109494 


3.67109494 


3.67109494 


3.67109 


0.1 


4.6288828 


4.6288828 


4.62888281 


4.62884 


0.3 


5.79657376 


5.79657363 


5.79657363 


5.79679 


0.5 


6.578402 


6.578402 


6.57840195 


6.57953 


0.7 


7.193266 


7.193265 


7.19326528 


7.19549 


1 


7.942405 


7.942404 


7.94240399 


7.94630 


2 


9.727325 


9.727322 


9.72732319 


9.73596 


50 


27.192660 


27.192638 


27.1926458 


27.2473 


1000 


73.419158 


73.419089 


73.419114 


73.5805 


8000 


146.745600 


146.745461 


146.745512 


147.0714 
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Table 4: Three-dimensional energies for states with k — 1 and I — 1,5, 10 
for V(q) = \ + ag 4 . Only £[4,4]'s are listed for different anharmonicities. 



a 


I = 1 


I = 5 


I = 10 


0.01 


4.76645813712 


9.289594583372 


15.233049583486 


0.1 


6.176138 


12.89579856 


22.309686916 


0.5 


8.93090 


19.3542918 


34.30436531 


1 


10.83313 


23.7006578 


42.25455311 


50 


37.4108 


83.258353 


149.9642236 


1000 


101.07403 


225.231013 


405.9901767 
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Table 5: Two-dimensional energies for states with k — 0, 1 and I = 0, 1, 5, 10 
for V(q) = *2+ aq A . Only £[4,4]'s are listed for different anharmonicities. 





a 


1 = 


/ = 1 


I = 5 


I = 10 





0.01 


1.0191783021 


2.056555600 


6.372257220 


12.0962676139707 




0.1 


1.150188 


2.4143403 


8.29606606 


16.976887733 




0.3 


1.33966 


2.895905 


10.53678440 


22.227347003 




0.5 


1.4760 


3.231453 


12.01658310 


25.611647809 




0.7 


1.5866 


3.499749 


13.17451055 


28.235574336 




1 


1.7242 


3.830324 


14.58077151 


31.403160969 




50 


5.512 


12.6399 


50.376652 


110.564887242 


1 


0.01 


3.129048426 


4.21691935465 


8.71236579303 


14.62528496652399 




0.1 


3.876642 


5.3954269 


12.01280444 


21.323943303 




0.3 


4.8105 


6.80426 


15.6275228 


28.28370297 




0.5 


5.4412 


7.74139 


17.9691075 


32.72664974 




0.7 


5.9389 


8.47655 


19.7879425 


36.15909940 




1 


6.5466 


9.3708 


21.9862477 


40.29318318 




50 


22.267 


32.237 


77.122811 


142.8926171 
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